Statistical Model of Evaporating Multicomponent Fuel Drops 

This model overcomes a deficiency of a prior statistical model. 
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An improved statistical model has 
been developed to describe the chemi- 
cal composition of an evaporating mul- 
ticomponent-liquid drop and of the 
mixture of gases surrounding the drop. 
The model is intended for use in com- 
putational simulations of the evapora- 
tion and combustion of sprayed liquid 
fuels, which are typically mixtures of as 
many as hundreds of different hydro- 
carbon compounds. Since an exact 
model providing a detailed account of 
all of the compounds would be compu- 
tationally intractable, the present statis- 
tical model is an approximation de- 
signed to afford results that are 
accurate enough to contribute to un- 
derstanding of the simulated physical 
and chemical phe- 
nomena, without im- 
posing an unduly 
large computational 
burden. 

As in any model of 
physical and chemical 
phenomena, some 
simplifying assump- 
tions are made: A 
drop is taken to be a 
sphere of radius R, 
wherein the liquid 
has constant density 
Pj. Evaporation of the 
liquid is assumed to 
occur under thermo- 
dynamic equilibrium. 

The gas surrounding 
the drop — a mixture 
comprising a carrier 
gas plus multicompo- 
nent vapor from pre- 
viously evaporated 
drops — is assumed 
to obey the perfect-gas equation of state. 
The gas is postulated to be quasi-steady 
with respect to the liquid, in the sense 
that the characteristic time of the gas is 
much shorter than that of the liquid. 
Consistent with what would be done in 
computations involving a very large 
number of drops, gradients within the 
drop are neglected; attention is paid 
only to volumetrically averaged proper- 
ties represented by the drop tempera- 
ture and the mass fractions of the chem- 
ical species in the drop. The focus on 
volumetrically averaged drop properties 
precludes consideration of phenomena 


associated with differences among diffu- 
sivities of different species. 

It is further assumed that the simu- 
lated phenomena occur at atmospheric 
pressure, where solubility of the carrier 
gas into the liquid is negligible, and that 
the far field conditions are quiescent. 
The model includes the applicable equa- 
tions for the conservation of mass, 
species, and energy. 

The statistical aspect of the model en- 
ters through invocation of the concept 
of continuous thermodynamics, accord- 
ing to which the chemical composition 
of a fuel is described probabilistically, 
by use of a probability distribution func- 
tion (PDF). This concept was summa- 
rized in two prior NASA Tech Briefs arti- 


cles: “Model of Mixing Layer With 
Multicomponent Evaporating Drops” 
(NPO-30505), Vol. 28, No. 3 (March 
2004), page 55, and “Simulations of 
Evaporating Multicomponent Fuel 
Drops” (NPO-30641), Vol. 29, No. 3 
(March 2005), page 72. However, the 
present statistical model differs from 
the model of the cited prior articles. 

In the prior model, the PDF is a sin- 
gle-peaked Gamma distribution, which 
is a function of the molar weight and of 
several parameters. However, the prior 
model does not generate a second peak 
at the low-molar-weight end of the 


drop-composition PDF that emerges 
when condensation of vapor onto drops 
occurs. 

The present model generates the 
needed second peak. The PDF in the 
present model is a superposition of two 
T-PDFs and, accordingly, is denoted a 
double-T-PDF. It is a function of the 
molar weight plus five parameters. 

Formally, the calculation prescribed 
by this model can be reduced to an in- 
verse mapping from (1) the First Five 
distribution moments, which can be 
calculated by use of the conservation 
equations, to (2) the five parameters of 
the double-T-PDF. The complexity of 
the inverse mapping, and the fact that 
the statistics of a discrete-chemical- 
species model may 
not be represented 
exactly by the dou- 
ble-T-PDF, make this 
calculation only ap- 
proximate. 

In a practical cal- 
culation, one makes 
a further approxima- 
tion by utilizing only 
the first four mo- 
ments and four pa- 
rameters plus a fifth 
parameter that is de- 
termined empirically. 
Despite these ap- 
proximations, exten- 
sive tests of the 
model on both diesel 
oil and gasoline show 
that the double-T- 
PDF accurately rep- 
resents the predic- 
tions of the 
discrete-chemical- 
species model (see figure). Moreover, 
the mean and variation of composition 
at the drop surface as computed by use 
of the double-T-PDF are in excellent 
agreement with those computed by use 
of the discrete-chemical-species model: 
this is an important result because these 
mean and variance determine the com- 
position of the gas mixture. 

This work was done by Kenneth Harstad, 
Patrick Le Clercq, and Josette Bellan of Cal- 
tech for NASA’s Jet Propulsion Laboratory. 
Further information is contained in a TSP 
(see page 1). 
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Double-r-PDF and Discrete-Model-PDF values were computed for a drop of diesel fuel evapo- 
rating at a gas temperature of 600 K. The results shown here are for the time at which a drop 
has evaporated to one-fifth of its initial mass. 
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